Tunneling between edge states in a quantum spin Hall system 
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We analyze a quantum spin Hall (QSH) device with a point contact connecting two of its edges. 
The contact supports a net spin tunneling current that can be probed experimentally via a two- 
terminal resistance measurement. We find that the low-bias tunneling current and the differential 
conductance exhibit scaling with voltage and temperature that depend nonlinearly on the strength 
of the electron-electron interaction. 
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A rapidly growing branch of condensed matter physics 
draws on the exploration of topologically nontrivial quan- 
tum states. Experimentally realized examples, which are 
by now well-understood, are given by the integer [l[ and 
fractional Q quantum Hall states. These states defy a 
classification in terms of the standard Ginzburg-Landau 
theory of symmetry breaking and a local order param- 
eter, but can instead be characterized by a topological 
quantity [l|, 0] • The importance of being able to identify 
a phase of quantum matter that does not fall under the 
Ginzburg-Landau paradigm has set off a search for other 
topologically nontrivial states, analogous to, but distinct 
from those connected to the quantum Hall effects. 

Some time ago, Kane and Mele — building on work by 
Haldane 0] — discussed the possibility of a new type of 
"topologically ordered" state of electrons in two dimen- 
sions: a quantum spin Hall ( QSH) insulator, proposed 
to be realized at low energies in a plane of graphene due 
to spin-orbit interactions [4j. Being a band insulator, a 
QSH insulator has a charge excitation gap in the bulk, 
but at its boundary there are gapless edge states with en- 
ergies inside the bulk gap. These states, which come in 
an odd number of Kramers' doublets, are "helical" (with 
clockwise/counterclockwise circling states carrying spin 
up/down, or vice versa, depending on the orientation of 
the effective electric field that enters the spin-orbit inter- 
action) and are responsible for the intrinsic spin Hall ef- 
fect that Murakami et al. had earlier predicted may occur 
in bulk insulators [Bj]. Time-reversal invariance implies 
that the energy levels of the counter-propagating edge 
states cross at particular points in the Brillouin zone. It 
follows that the spectrum of a QSH insulator cannot be 
continuously deformed into that of an ordinary band in- 
sulator, which has zero (or equivalcntly, an even number 
of) Kramers' doublets. In this exact sense, a QSH insula- 
tor realizes a topologically nontrivial state of matter Q. 
In subsequent and independent work, the QSH insulator 
state was proposed to occur also in strained semiconduc- 
tors and in HgTe quantum wells with an "inverted" 
electronic gap 0]. An experiment carried out on quasi- 
two-dimensional HgTe quantum wells grown by molec- 
ular beam epitaxy and sandwiched between (Hg,Cd)Te 
barriers has revealed data consistent with helical edge 
state transport, suggesting the first observation of the 



QSH effect The possibility of dissipationless trans- 
port of spin currents along the edges of a QSH insulator 
is a tantalizing prospect for future spintronics applica- 
tions [Tol | - T° make progress, however, a more complete 
picture of the physics is required. 

An important issue is to understand the behavior of 
edge currents in the presence of a tunneling junction con- 
necting two opposite edges of a QSH bar (FIG. 1). When 
the bar is connected to a battery, a net spin current can 
tunnel through the junction, and one would like to know 
how the electron-electron interaction influences its con- 
ductance. This is the problem we shall address here. 
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FIG. 1: (Color online) Geometry of the QSH point contact 
device studied in this paper. The full (dotted) lines represent 
helical edge states in equilibrium with the left (right) contact. 

We consider the simplest situation with a single 
Kramers' doublet of helical edge states, applicable to a 
tunneling experiment on a HgTe quantum well 0| • In the 
absence of electron interactions this case can formally be 
thought of as resulting from a superposition of two inte- 
ger quantum Hall systems with the up- and down-spins 
of the electrons being subject to opposite effective mag- 
netic fields. This emulates the spin-orbit interaction that 
is built-in in the k-p Hamiltonian that defines the elec- 
tron dynamics close to the Fermi level of the quantum 
well [9J. The bar is connected to a battery with left (L) 
and right (R) contacts as in FIG. 1. Applying a gate 
voltage V g perpendicular to the upper and lower edge 
of the bar at x = will bring the edges close to each 
other, forming a point contact at which electrons may 
tunnel from one edge to the other. With the gate volt- 
age turned off there is no tunneling present, assuming 
the edges to be well separated. For the case illustrated 
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in FIG. 1, electrons originating in the L [R] contact of 
the battery carry current to the right [left] , with spin-up 
[spin-down] on the upper edge and spin-down [spin-up] 
on the lower edge. The right- [left-] moving electrons are 
in equilibrium with the left [right] contact and have a 
Fermi energy equal to the electrochemical potential 
[Hn] of that contact. Note that counterpropagating elec- 
trons do not equilibrate when injected at different chem- 
ical potentials since any scattering off impurities or de- 
fects conserves spin, thus making impossible transfer of 
electrons from one type of edge state to the other. If the 
driving voltage V= ([J,i,—fiR.)/e > 0, a net charge current 
flows from left to right on each edge, accompanied by a 
spin current carrying spin-up [spin-down] on the upper 
[lower] edge. Neglecting electron interactions, the ratio of 
the drain-source charge [spin] edge current to the driving 
voltage is the Hall [spin Hall] conductance e 2 /h [e/4ir] 
[111 ]. When the gate voltage V g is turned on, more of the 
right-moving electrons tunnel through the point contact 
at a finite driving voltage V > 0, leading to a depletion of 
the source-drain current. While there is no net tunneling 
of charge between the edges, the point contact supports 
a net inter-edge spin tunneling current (cf. FIG. 1). 

The picture becomes more complex when allowing for 
electron-electron interactions at the edges. Away from 
half-filling of the one-dimensional (ID) band of edge 
states, time-reversal invariance constrains the possible 
scattering processes at an edge to dispersive (d) and for- 
ward (/) scattering (THH. In the vicinity of the Fermi 
points the corresponding interactions are given by 

H d =9d J dx (i'\ l ^ m ^l l ^ Ll + ipl^L^^-ipR^ (1) 

and 

H f=9f / dx ^LVw^LVw (2) 

a = R,L 
" = 1.l 

Here ipRi and ipL\ are ID fields that annihilate an elec- 
tron in a clockwise propagating helical state on the up- 
per and lower edge, respectively. Similarly, tp^i and V'-RI 
are fields that correspond to a counterclockwise propa- 
gating state on the upper and lower edge, respectively. 
It is here important to emphasize that the presence of 
the dispersive scattering channel, controlled by (P), is 
a fundamental difference between the edge physics of a 
QSH insulator and a system exhibiting the integer quan- 
tum Hall effect (IQHE). As a result, the QSH insulator 
may show interaction effects which arc suppressed in the 
IQHE. Adding a linearized kinetic term 

H = -iv F ^ / dx {^ R(T d x ip Ra - V'l^^Vio-) i (3) 
we can bosonize H = H + + Hf, and obtain 

H = \ E [ dx {U d ^f + K{d x eA , (4) 

1=1,2"' ^ ' 



^ K = y lSppS and v = vWf)Mf) 2 , 

vf being the Fermi velocity. The indices 1 and 2 label 
the upper and lower edge, respectively, with 4>i = if> R ^ + 

4>li,4>2 = 4>li:+4>ri^i = <Art- 0i4 > and 6» 2 = M-<t>Ri> 

where 4> aa and 8 aa define chiral boson fields and their 
duals within the standard bosonization scheme [3] ■ Note 
that in contrast to an ordinary spinful Luttinger liquid 
which exhibits spin-charge separation, the boson fields in 
(U]) contain both charge and spin. While hclicity makes 
spin a redundant quantum number on a single edge, it 
is important to include it when two edges are connected 
via a point contact. The tunneling through the contact, 
with amplitude u, is governed by the operator 

Ht = u (VltV'fiT + iPr^LI + i>Rii>Li + ^li^Rl) i ( 5 ) 

defined at x — 0. It can similarly be bosonized: 
2u 

H t = — sin[ViOi + h)} coslv^i + 6 2 )\. (6) 

Given the bosonized theory, eqs. (jU) and © , we may now 
use standard pcrturbative RG arguments to uncover the 
effect of electron interactions on the tunneling. 

As a first step, we integrate out the bosonic fields in 
the partition function of the system except at a: = 0, 
thus obtaining a theory defined only at the location of 
the point contact [l5|. With A an energy cutoff, and 
t = it Euclidean time, this gives 

Z ~ / J] exp(~S - S t ), (7) 

i=l,2 

where 
and 

S t = -^J dr sin [V^ (Mr) + Mr))] 

xcos[V^(0i(t)+0 2 (t))]. (9) 

Next, the localized fields are split into slow (s) and fast 

(/) modes, Mr) = J% <t^ T M) and ^(r) = 

f\/b<\oj\<AlTT e ~~ lUJT< t ) (. UJ )> with b > 1 a scale factor, and 
with a similar definition of 6i S and 9if (i = 1,2). A 
cumulant expansion in u then gives an expression for the 
low-energy effective action, call it S e ff. To 0(u 2 ), 

e -s eff [t,] =e -s 3 [0 s ] e <s t > / -i((s t 2 ) f -(s t >J.)+..^ (1Q) 

Here S s [<p s ] is the slowly fluctuating part of S, while 
(. . .)j is an average taken over the fast modes. The cal- 
culation of (St) f and the second-order cumulant (Sf) , — 
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(S t ) 2 f is here somewhat cumbersome, but is facilitated by 
the presence of the time-reversal symmetry. We find that 



(S t ) f = -b-^ K +WS t i<j> ls ,<l>2 S ,6u 



fill 



with St[4>is, 4>2s, 9isi &2s] as m © but with the slow fields 
replacing the original ones. As for the second-order term, 

(S 2 ) f - (S t ) 2 f = ^ J dr (V e cos {29 1 + 20 2 )] 

+ l/ cos[V^ (20i+20 2 )] +...), (12) 

where V e = b 1 ~ 2 / K -b l ~ K - 1 / K , V = ^-K-XjK _ h x-iK ^ 
and where . . . indicate higher-order terms that do not 
influence the renormalization to this order in u. The 
first-order RG equation for u, 



du ( 1. 1 

= u 1 (A H ) , 

dm 6 \ 2 V K n 



(13) 



is obtained from (fT0| and (fTT|) and reveals that the scal- 
ing dimension Ak of the tunneling operator Ht in j6|) is 
Aa' = \{K + 1/A). As for the second-order equations, 
these are extracted from (TIT)]) and (Tl"2"|) , and read 

dVe__u 2 f M 2 
d\nb 



4 f(l - l) e ( 1 - 2 /^) i" b -(i-2A K )e^- 2A ^ ln b ) 
7r y A y 



^ =4 f(l-2A K )e( 1 - 2A - )lnb -(l-2A) e 
a In o 7r V 



l-2/f)lnfc 



These equations imply that, to second order in it, H t 
renormalizes to zero for all values of K in the interval 
1/2 < K < 2. This includes the experimentally rele- 
vant regime for a HgTe quantum well: A rough estimate 
of A for this case, based on the approximate relation 
Km(l + U/{2E F ))~ 1 ' 2 0, yields that 0.8 < A < 0.9, 
using that Ep w h/m*r 2 and U ~ e 2 /er s , where e is 
the electron charge, e w 20eo is the dielectric constant, 
m*«0.02m e (with m e the electron mass) [13], and where 
r s is the effective Bohr radius for electron densities n e in 
the interval 0.5 x 10 11 cm~ 2 < n e < 3.5 x 10 11 cm -2 (at 
which the experiment in Ref. was carried out). 

It is interesting to compare the "weak-tunneling" fixed 
point found here to the situation for the quantum Hall 
effects, where for the IQHE the tunneling between edge 
states is marginal, while for the fractional quantum Hall 
effect the tunneling renormalizes to large values for all 
filling fractions In contrast, as seen from the second 
order RG equations above, a "strong-tunneling" QSH 
fixed point appears only for K < 1/2 (or for the unphys- 
ical region K > 2 with attractive electron interaction). 

Turning to the tunneling current in the presence of a 
driving voltage V = — [iRj/e, we shall focus on the 
case of a low bias, allowing us to use a linear response 
formalism The current I c (t) that we shall calculate 
is the sum of the charge tunneling currents between edge 
states with the same helicity, related to the total spin 



tunneling current I s (t) by I c (t) = (2e/H)I s (i). Since I c (t) 
is equal to the depletion of the charged source-to-drain 
current in the presence of the point contact (cf. FIG. 1), 
it follows that the spin current I s (t) can be detected ex- 
perimentally by a two-terminal resistance measurement. 

With V > 0, the current I c (t) can be expressed as the 
rate of change of the number of electrons in equilibrium 
with the left contact of the battery (see FIG. 1), I c (t) = 
— e(iVi(t)). The number operator Nl = a(i/>j ? |'0i?.i + 
ipR^ipRi), where a is a lattice constant, has the property 
that Nl = i[H+H t ,N L ] = i[H t ,N L ]- This implies that 

I c (t) =eu 2 J dt'®{t-t')(e le tt' dt " v{ - t "> (\A{t),A\t')]) 
_ e -ie /*' dt"v(t") ^[ jA t( t ) i A (t')]) ) , (14) 



where A = a(ip L *ipitf + ip L ,ipui)- Introducing the re- 
tarded Green's function G ret (t) = -iO(t) ([A(t), A^(0)]) 
and its transform G re t(— eV) = J dte~ zeVt G re t{t), it fol- 
lows that for constant V the integral in eq. (fT4|) can be 
written as — 2Im[Gr re t(— eV)]. The correlation functions 
G+(t) = (A(t)A*(0)) and G_(t) = (A*(0)A(t)) are eas- 
ily calculated in the bosonized theory, and one finds that 



G±{t) - \ \-v(t±iS) 



2Aj 



(15) 



where 5 is a short-time cutoff. Collecting the results, 

\2A K 



I r = 2eu 2 ^\ , (eV) 2A > 



T(2A K ) - ' (16) 

which tells us how the dc tunneling current scales with V 
in the limit V — > 0, and also how its amplitude depends 
on the parameter A that encodes the electron interac- 
tion. To account for the full A-dependence in eq. p^|) one 
uses the parameterizations of A and v after eq. (|4|), with 
9d w 4 5/ 0. To 0(g f /v F ), v w v F (5 + 3A)/(3 + 5A). 

In order to extract the finite-temperature tunneling 
conductance G, we perform a conformal transforma- 
tion of the correlation functions in (|15[) . first going 
to Euclidean time r = it, and then taking vt — > 
(v/3/2tt) arctan(27rr//3w), with 0=1/T. It follows that 



I c = -2eu 2 (a/vy* K (2irT) 



\2A f 



x Im 



B(A K + ieV/2irT, A K - ieV/2irT)} 
sin(7r(A jr - ieV/2irT)) 



(17) 

cos(ttAa-) 

where B is the Euler beta function. In FIG. 2, the current 
is plotted for a few different values of A and T. We have 
here taken a w 1 nm and vf ~ 6 x 10 6 m/s SEEl, 

and 

put u = O.lvp/a. From (|17|) we obtain the scaling of the 
zero-bias conductance G with temperature T, 



r _ dl c 
G= dV 



cx T 



2A K -2 



(18) 
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FIG. 2: (Color online) The two graphs show the charge tun- 
neling current I c as a function of the applied voltage Vfor dif- 
ferent values of K and T . (The spin tunneling current I a that 
transfers spin between the edges is given by I s — (h/2e)I c .) 

It is also interesting to explore the tunneling current 
for an ac voltage of the form V(t) = Vq + V\ sin(f^). 
Inserting V(t) into eq. (fT4"|) and following Ref. [2CJ, we 
find that the dc component / c ,o of the current, defined 
as the time average of I c (t), can be expressed as / Cj o = 
2eu 2 {a/v) 2AK J2n a n(eVi/n)(eV + 7ift) 2A ^-\ where 

1 f27r f 27r ^^/gin(t'-i) gi^L ( cos t'-cos t) _ 



o»(e^i/fi) 



(2tt) 



Jo Jo 



In FIG. 3 we have plotted the dependence of I c ,o on Vb 
for some different values of K and V\. As seen from the 
figures, / c .o decreases with increasing electron-electron 
interaction (i.e. with decreasing values of K). 
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FIG. 3: (Color online) The dc component of the charge tun- 
neling current 7 Cj o as a function of Vb for different values of 
K and Vi. (The dc component of the accompanying spin 
tunneling current I s ,o is given by I 3t o — (fi/2e)/ c .o.) 

To summarize, we have found that a point contact con- 
necting two edges of a QSH bar supports a spin tun- 
neling current / cx v 2Ak ~ 1 at small voltages V, with 
a zero-bias conductance G cx T 2Al< ~ 2 for all values of 
A K = (K + l/K)/2 with 1/2 < K < 2, where K en- 
codes the strength of the electron interaction. This spin 
current can be probed experimentally via a two-terminal 
resistance measurement. The interval 1/2 < K < 2 con- 
tains the if-values applicable to a HgTc quantum well in 
the QSH regime When K < 1/2, the tunneling ampli- 
tude scales to large values, effectively severing the edges, 
analogous to what happens in a fractional quantum Hall 
system. Given that a QSH device can be manufactured 
which allows K to pass through the value of 1/2, this 
would open for the possibility to experimentally study 



the transition between strong and weak tunneling in a 
topologically nontrivial phase of matter. 
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Note added: Upon completion of this work we found a 
preprint by Hou et al. (arXiv:0808.1723vl, published in 
Ref. [HI), on QSH edge states in a four- terminal corner 
junction geometry. For weak and intermediate electron 
interactions these authors find a weak-tunneling fixed 
point, similar to our result for a two-terminal device. 
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